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We experimentally investigate the magnetic field B dependence of the critical current I c and the 
transition temperature T c , i.e. the I C (B) —T C (B) phase boundary, of superconducting niobium thin 
films patterned with periodic and quasiperiodic antidot arrays on the submicron scale. For this 
purpose we monitor current- volt age characteristics at different values of B and T. We investigate 
samples with antidots positioned at the vertices of two different tilings with quasiperiodic symmetry, 
namely the Shield Tiling and the Tuebingen Triangle Tiling. For reference we investigate a sample 
with a triangular antidot lattice. We find modulations of the phase boundary for both quasiperiodic 
tilings, which were predicted by numerical simulations but not observed in experiments yet. The 
particularity of these commensurability effects is that they correspond to excess flux densities, 
which are slightly higher than the matching flux. The observed matching effects can be explained 
by quasiperiodic caging of interstitial vortices and/or the formation of symmetry induced giant 
vortices. 

PACS numbers: 74.25.Uv, 74.25.Sv 74.25.Dw, 75.50.Kj 



I. INTRODUCTION 

The topology of superconducting thin films has a 
strong impact on the shape of the superconductor- 
normalconductor phase boundary in an external mag- 
netic field B. Basically the reason for this connection is 
the fluxoid quantization condition, requiring that along 
each closed path inside the superconductor the gauge in- 
variant phase is an integer multiple of 2ir. 

As pointed out first by Little and Parks [ ] this quan- 
tization condition leads to an oscillation of the transition 
temperature T c of a superconducting cylinder in a mag- 
netic field. In a field parallel to the cylinder symmetry 
axis, the periodicity of the T c (£?) -oscillations is one flux 
quantum <£> = h/2e = 2.07 x 10 -15 Tm 2 per cross sec- 
tional area of the cylinder. These Little-Parks oscillations 
were also found in two-dimensional (2D) superconducting 
wire networks with periodic or quasiperiodic structure 
[2-4] , where the flux quantization has not only to be ful- 
filled for each single loop of the network but also for all 
paths around multiple loops. In these samples, commen- 
surable states between the applied flux and the network 
are also observable at fractional (periodic networks) or 
even irrational (quasiperiodic networks) multiples of one 
flux quantum per single loop. 

Another famous consequence of the fluxoid quantiza- 
tion is the formation of Abrikosov vortices, when a type- 
II superconducting thin film is exposed to a perpendic- 
ular magnetic field. In an ideal type-II superconductor, 
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these vortices with normal conducting cores are repul- 
sively interacting and form a highly ordered state, the 
hexagonal Abrikosov vortex lattice [ ] . As the formation 
of a normal conducting core costs condensation energy, 
any defect in a real superconductor, where the supercon- 
ducting phase is weakened, can be viewed as an energy 
minimum for an Abrikosov vortex. For the case of an 
ensemble of vortices interacting with an ensemble of de- 
fects, the repulsive interaction of the vortices plays an 
important role and the resulting spatial vortex configu- 
ration is determined by an interplay of the energy gain 
due to vortex pinning at defects and the energy cost due 
to the elastic vortex lattice deformation. 

With modern lithography techniques it is possible to 
create lattices of artificial pinning sites, e.g., antidots 
[6, 7], magnetic dots [ , ] or carbon nanotubes [10] with 
almost arbitrary topology, on the characteristic length 
scales of the superconductor (penetration depth A and co- 
herence length £). Such pinning potential landscapes can 
be varied in terms of size, geometry and amplitude dur- 
ing the patterning process. In addition, the strength and 
range of the vortex- vortex as well as the vortex-defect in- 
teraction can be tuned during the measurements by vary- 
ing the sample temperature T and hence A and £. Finally, 
it is possible to easily control the vortex density via the 
applied magnetic field and to create constant and alter- 
nating forces on the vortices by applying external cur- 
rents. Hence, Abrikosov vortices in a micropatterned su- 
perconducting film constitute a highly designable model 
system for interacting particles in a 2D potential. 

In such systems, dynamic effects like phase locking 
phenomena [ 1-13] and ratchet effects [14, 15] have been 
intensely investigated, as well as static effects like the for- 
mation of vortex quasicrystals [16-20] in particular for 



2 



artificial defects arranged in a Penrose pattern. An ex- 
perimentally easily accessible quantity to gather integral 
information on the static interaction between a vortex 
lattice and a pinning array is the magnetic field depen- 
dent critical depinning current I C (B), where e.g. com- 
mensurabilities between pinning sites and vortex lattices 
can be seen as pronounced maxima [6-9, 16-20]. 

The implementation of artificial defects as pinning sites 
for Abrikosov vortices has furthermore relevance for a va- 
riety of applications, as in many cases the dissipative 
vortex motion leads to a reduction or a limitation of 
the performance of superconducting devices and hence 
is desired to be suppressed by pinning sites. Approaches 
for improving the performance of superconducting mi- 
croelectronic devices have found strategically positioned 
antidots feasible for the reduction of low-frequency flux 
noise in quantum interference devices [ ] and for the re- 
duction of vortex associated losses in coplanar microwave 
resonators [22, 23]. Quasiperiodic arrays might be par- 
ticularly suitable for some of these applications, as they 
have many built-in periodicities, which moreover are not 
only found for integer or rational but also for irrational 
multiples of the matching field, i.e. the field for which 
the vortex density n v equals the density of pinning sites 
n p . For the Penrose tiling for instance, this leads to a sig- 
nificant broadening of the I C (B) peaks, as compared to 
the triangular defect lattice [ - ] and hence to a more 
homogeneous I C (B) dependence. 

In this work, we experimentally investigate the I C (B) — 
T C (B) phase boundary of superconducting niobium thin 
films with three different arrangements of antidots - two 
kinds of quasiperiodic arrays, namely the Shield Tiling 
(ST) [ ] and the Tuebingen Triangle Tiling (TTT) [25] 
- and for reference a triangular lattice. We perform our 
measurements close to T c , where Little-Parks oscillations 
and collective Abrikosov vortex pinning are not strictly 
separable effects. There might also be a smooth transi- 
tion with decreasing T from the wire network limit to a 
film with holes [26, 27]. However, we are not going to fo- 
cus on such a possible transition in this paper, but rather 
investigate commensurable states between an ensemble of 
quantized fluxoids and a given microtopology. Hence, we 
monitor the phase boundary holistically. We find and 
discuss signatures for commensurabilities at nonmatch- 
ing excess flux densities, which have been predicted by 
theory [28] and which might be attributed to a caging 
effect and the formation of giant vortices. 

The paper is organized as follows. After the introduc- 
tory Sec. I, we introduce the sample fabrication in Sec. II 
and the method for their characterization in Sec III. In 
Sec. IV we present and discuss the experimental results 
on the phase boundary of niobium thin films with anti- 
dot arrays arranged in a Shield Tiling and a Tuebingen 
Triangle Tiling, respectively. We also discuss experimen- 
tal results on a sample with a hexagonal antidot array 
for reference. Section V contains the conclusion of this 
work. 



II. SAMPLE FABRICATION 

The experiments were carried out on cross-shaped 
structures patterned in d = 60 nm thick dc magnetron 
sputtered Nb films. Figure 1 (a) shows a sketch of 
a 800 x 800 fim 2 sample (white areas are Nb) with a 
100 x 100 fim 2 large center area (bridge), which is pat- 
terned with different arrangements of circular antidots. 
Independent of the specific arrangement and array sym- 
metry, the antidot density in the center area of these 
bridges is n p « 2 /mi -2 corresponding to a total num- 
ber of antidots N p ~ 2 • 10 4 and matching field Bi = 
rip&o ~ 4mT. The patterning of the structures including 
the antidots was performed by e-beam lithography and 
subsequent SF6 reactive ion etching. 
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FIG. 1: (Color online) (a) Layout of a 800 x 800 /nn 2 chip 
with a cross shaped Nb structure for a four-probe current (/)- 
voltage (V) characterization of a 100 x 100 /mi 2 large area pat- 
terned with antidots (center square), (b)-(d) Sketches of the 
tilings and atomic force microscopy (AFM) images of the an- 
tidot arrays showing (b) the triangular/hexagonal lattice, (c) 
the Shield Tiling and (d) the Tuebingen Triangle Tiling; the 
antidot density of all three arrays is n p « 2 /im" 2 and the 
diameter of the circular antidots is D w 300 nm. 

Each bridge contains either a triangular/hexagonal ar- 
ray of antidots or quasiperiodically arranged antidots at 
the vertices of the Shield Tiling or the Tuebingen Tri- 
angle Tiling. The three different patterns are shown in 
combination with a sketch of the tiling and an atomic 
force microscopy (AFM) image of the corresponding an- 
tidot array in Figs. 1 (b), (c) and (d). A more detailed 
description and discussion of the quasiperiodic tilings is 
given together with the results in Sec. IV. The antidots 
of the samples shown here and used for the presented 
experimental results have a diameter D « 300 nm. 

For a rough estimate of the magnetic penetration 
depth A and the coherence length £ in our samples, we 
determined the resistance slightly above T c as R n « 



FIG. 2: (Color online) Data for Nb structure with triangular antidot lattice at variable temperature T and magnetic flux 
density B (normalized to the matching field B\ « 4mT). (a) Critical current I c for a voltage criterion V c — 5 pN\ (b) Resistance 
R measured with a bias current / = 50 fiA. 



0.6 (cf. Sec. IV) and calculated the resistivity p = 
R n A/L = 3/i(]cm with the center area cross section 
A — 100 pm x 60 nm and length L = 100 pm. From 
simulations we know that due to our cross-shaped ge- 
ometry this resistivity has to be corrected by about a 
factor of two, as the bias current / is not confined ex- 
actly to the center area. So for further calculations we 
use the corrected resistivity p = 2 • p = 6 pQcm. Using 
pi = 3.72 x 10~ 6 pQcm 2 [ ] we determine the electron 
mean free path / « 6nm. As we have neglected the 
antidots in this calculation, the value for p is probably 
somewhat over-estimated and as a consequence the value 
for I is somewhat under-estimated. With the BCS coher- 
ence length for niobium of £o = 39 nm and the "dirty 
limit" expression [ ] f(T) = 0.855>/Jfo/ ' y/l - T/T c we 
find £(T) = 13 nm/ ^/l — T/T c . Analogously, with the 
London pe netration dept h A^ = 38 n m at T = and 
A(T) = VV1.33ZA L / V 2 (! - T/Tc) we find A(T) = 
59 nm/^! - T/T c . As for all our data T/T c > 0.98 and 
hence A(T) > 300 nm^ d = 60 nm, we have to consider 
the thin film penetration depth A = 2X 2 /d. Again, £(T) 
and A(T) are somewhat under- and over-estimated, re- 
spectively. 



III. SAMPLE CHARACTERIZATION 

The samples were mounted in a low-temperature setup 
with high temperature variability 4.2 K< T < 100 K and 



stability ST < lmK [ ]. We record current- volt age 
characteristics (IVCs) over the whole 2D phase space of 
perpendicular magnetic field between the critical fields 
±B C 2 and temperatures close to the transition tempera- 
ture T c . We keep the stepwidth between different tem- 
peratures AT and different values of magnetic field A£> 
constant; so for each value of T we have IVCs for all val- 
ues of B and vice versa. This procedure enables us to 
extract all information as I C (B), T C (B) or R(B)-cuives 
for arbitrary and freely selectable voltage criteria V c , re- 
sistance criteria R c and bias currents I. 

Figure 2 shows for a sample with a triangular antidot 
array (a) the critical current I C (T,B) and (b) the re- 
sistance R(T,B). From Fig. 2 (a), we can extract single 
vertical and horizontal slices, which correspond to I C (B)- 
curves (dark lines) for constant T and to T c (£>)- curves 
(bright lines) for constant J, respectively. Similarly, from 
Fig. 2 (b) we get R(B)-curves for constant T (dark lines) 
and T c (£?) -curves for constant V (bright lines) by taking 
vertical and horizontal slices, respectively. 

The cross shape of our Nb structures is expected to 
induce an inhomogeneous current density distribution 
across the center area. This point has been particularly 
discussed for a ratchet experiment [32-34] , where the di- 
rection of the local driving force relative to the intrinsic 
symmetry axes of the pinning lattice is of great impor- 
tance. Also in periodic and quasiperiodic potential land- 
scapes with circular antidots, the vortex dynamics might 
depend on the direction of the local transport current 
density [35-38]. However, we do not find our experi- 
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mental results to differ significantly for the two possible 
perpendicular directions of the bias current. It is most 
likely anyway that the depinning current density is first 
reached directly between the voltage pad contact points, 
where the mean current density reaches its maximum due 
to the smallest cross section area. So we think that the 
essential results of this work are not influenced by the 
bridge geometry. 

Due to the used voltage amplifier, we had to face a 
voltage noise of several 1CT 7 V. To improve the signal to 
noise ratio and the visibility of small matching features, 
the results presented in this paper are not single-shot raw 
data but have been averaged and smoothed for their final 
interpretation and presentation. First, we measured five 
IVCs at each B — T phase space point, each consisting of 
1000 data points, and averaged them for the further pro- 
cessing. The averaged IVCs were slightly smoothed by a 
non weighted five-point adjacent- aver aging before we ex- 
tracted voltages for pre-defined currents or vice versa. At 
the end we once more performed a non weighted three- 
point adjacent- averaging on the resulting I C {B), T C (B) 
and R(B) slices. 



IV. EXPERIMENTAL RESULTS 
A. The Triangular Lattice 

Before we discuss our results on quasiperiodic antidot 
arrays, we present the results on the sample with the tri- 
angular lattice. The data in Fig. 2 already reveal clear 
and strong features, which indicate commensur abilities 
between the flux line and antidot density. We find pro- 
nounced and narrow ridges in (a) and canyons in (b) at 
B — B\ (corresponding to one vortex per antidot). We 
find quite similar but somewhat weaker structures for flux 
densities corresponding to two and three vortices per an- 
tidot (B/Bi = 2 and BjB\ =3). Our findings show that 
our samples and our data analysis are suitable to reliably 
detect not only strong commensurability features at in- 
teger matching fields, but also weaker ones at fractional 
matching flux densities. These weaker commensur abili- 
ties at fractionals of the matching flux reveal themselves 
by characteristic fine-modulations of the phase boundary 
at flux densities between B = and B\. 

Figure 3 (a) shows several I c ( B) -slices. Marked with 
arrows between B = and there are four small sub- 
peaks visible, which we attribute to the fractional fill- 
ing factors 1/4, 1/3, 2/3 and 3/4 (dashed lines). Corre- 
sponding vortex configurations and in addition the one 
for B = B\ are depicted schematically in Fig. 3 (b). The 
observed fractional matching peaks for < B < B\ and 
the corresponding vortex configurations have been found 
in theoretical as well as some experimental studies on tri- 
angular pinning lattices before, cf. e.g. [39, 40]. We find 
very similar peaks for flux densities between B = B\ and 
B = 2B u i.e. for B/B 1 = 5/4, 4/3, 5/3 and 7/4, indicat- 
ing that the same configurations are stable for two- vortex 



occupations, that is when all antidots are already occu- 
pied by one vortex. 



B. The Shield Tiling 

The Shield Tiling [cf. Fig. 1 (c)] consists of three differ- 
ent tiles: a square, an equilateral triangle and a deformed 
hexagon with alternating interior angles of ir/2 and 57r/6. 
This is a very interesting collection of tiles, as two of 
them locally resemble the energetically most favorable 
configurations of a periodic Abrikosov vortex lattice: the 
triangular and square array. Each hexagon as the third 
basic shape leaves a rather large antidot-free area, what 

- regarding the whole lattice - leads to quasiperiodically 
distributed, individual cages or reservoirs for interstitial 
vortices as also described in Ref. [28]. The pattern was 
created by starting with six hexagonally arranged trian- 
gles followed by an iterative application of inflation rules 
[24]. 

In a previous theoretical study on different quasiperi- 
odic pinning arrays, we have found that the critical cur- 
rent density of a superconductor with antidots arranged 
in the Shield Tiling shows two pronounced I C (B) max- 
ima. One of them corresponded to the first matching field 
B = B\. The second maximum was found for an excess 
flux density of B « 1.18£?i. According to simulations, 
this corresponds to the situation, where all antidots are 
occupied by one vortex, and within each hexagon there 
is an additional interstitial vortex [ 28] . These interstitial 
vortices are caged, i.e. pinned by the repulsive interac- 
tion with the surrounding vortices, which are pinned at 
antidots. 

In the sample investigated here, we expect the addi- 
tional maximum at B « 1.16£?i, which is at a slightly 
lower value than in Ref. [ ]. This difference is due 
to the fact that the maximum is expected at B = 
(1 + Nh/N P )B\, where Nh and N p is the total number 
of hexagons and antidots, respectively. The ratio Nh/N p 
converges for large numbers of antidots. In the numerical 
approach in Ref. [28] we have only used N p ~ 500 pinning 
sites with Nh/N p « 0.18. In our experimentally used ar- 
ray with N p « 20000 antidots we have N h /N p « 0.16. 

Fig. 4 (a) shows 7 c (I?)-curves of the Shield Tiling sam- 
ple at several different temperatures. At first sight, there 
seems to be no additional maximum around B = l.WBi 
but only a very strong and single maximum of the critical 
current at B = B\. However, when we cool the sample to 
lower temperatures a double-peak structure appears with 
one peak at B = 2B\ and one at a slightly higher field 
value B « 2.16£?i. This additional peak - or more pre- 
cisely speaking its height relative to the one at B = 2B\ 

- has a strong dependence on temperature. It is almost 
absent at T = 8.413 K but then grows and becomes the 
dominating peak at T = 8.393 K (and probably below). 
A very similar relation and development between two 
peaks with temperature has been reported for caged vor- 
tices in periodic pinning arrays before [41] and might be 
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FIG. 3: (Color online) (a) Critical current I c vs normalized magnetic flux density B / B\ of a Nb structure with a triangular 
antidot lattice at five different temperatures T; dashed vertical lines indicate the fractional flux densities B / B\ — 1/ '4, 1/3, 2/3 
and 3/4 and arrows point to the corresponding local maxima, (b) Sketches of (possible) vortex configurations corresponding 
to local maxima of I c at Bj B\ — 1/4, 1/3, 2/3, 3/4 and 1. 



attributed to an increasing interaction strength between 
the vortices due to a decreasing penetration depth. 

The predicted matching at an excess vortex density 
peak seems to appear not at B = 1.1 6Bi but at B = 
2.165i. It is thus very likely that singly occupied anti- 
dots are still too attractive for vortices entering above 
B = Bi, hence, they first occupy the antidots twice 
before the centers of the hexagons become energetically 
more favorable. We show sketches of the proposed vor- 
tex configurations in Fig. 4 (b). Of course, the situation 
depicted here is somewhat idealized. In simulations we 
find similar situations, but we also find a certain degree 
of disorder in the vortex lattice in any case. For exam- 
ple, for B = 2Bi it might well be that not all six antidots 
of the hexagons are doubly occupied. It is rather likely 
that one or two of them only host one vortex, while there 
are also interstitials present. The simulations we refer to 
here, are not shown, as they were not performed exactly 
with the experimental parameters. However, they served 
as an inspiration for the proposed vortex configurations. 

The discussed observations of commensurabilities in 
the Shield tiling are closely related to widely discussed 
caging effects in different kinds of periodic pinning arrays, 



including randomly or periodically diluted triangular lat- 
tices [41-45]. These studies have revealed, that matching 
peaks can be shifted to higher field values than those, for 
which vortex and antidot density are equal. Whenever 
there are areas with missing antidots, these areas fill up 
with interstitials, when the antidots are saturated and 
these configurations with interstitials can be more stable 
than the configurations without interstitials. 



C. The Tuebingen Triangle Tiling 

The second quasiperiodic arrangement - the Tuebin- 
gen Triangle Tiling - is depicted in Fig. 1 (d). It consists 
of two different isosceles triangles, one acute-angled and 
one obtuse angled. This tiling can also be constructed 
by an iterative application of inflation rules [ ] after a 
starting configuration of ten acute-angled triangles. In 
contrast to the triangular and the Shield tiling, the edge 
lengths are not equal for all basic shapes here, but we 
find two different lengths, whose ratio is the golden mean 
r w 1.618. This can lead to a locally reduced pinning 
site density, when ten triangles are composed such that 



6 




BIB X 



(b) • Vortex • Double-Vortex 




B =B X B = 2B X B = 2A6B l 



FIG. 4: (Color online) (a) Critical current I c vs normalized magnetic flux density B / B\ of a Nb structure with a Shield Tiling 
lattice of antidots at five different temperatures T; dashed vertical lines indicate the fractional flux density B/Bi = 2.16 and 
arrows point to the corresponding local maxima, (b) Sketches of (possible) vortex configurations corresponding to maxima of 
the critical current at B / B\ — 1,2 and 2.16. 



10 antidots form a ring, for example in the left upper 
and lower corner of Fig. 1 (d). Below, these full rings 
with 10 antidots surrounding one as a circle are called 
10/10-circles. 

The 10/10-circles have a high antidot line density on 
the circle and a reduced density in the interior and may 
again be viewed as cages or reservoirs for interstitial vor- 
tices. In contrast to the cages in the Shield Tiling we 
here find the situation that there actually is one anti- 
dot in the center of each 10/10-circle. Analogously to 
the Shield Tiling discussion we have to calculate the 
ratio of the number of 10/10-circles AT 10 over the to- 
tal number of antidots N p in order to assign the ap- 
pearance of local maxima in I C (B) to vortex configura- 
tions related to the circles. This ratio is found to be 
Niq/N p = 1/t 6 ~ 0.056 for N p — >• oo, as after each itera- 
tion the total number of antidots increases by a factor of 
t 2 , i.e. = r 2 A^ _1 for the ith iteration. After three it- 
erations each antidot has transformed into a 10/10-circle, 
i.e. N{ = N*- 3 = (l/r 6 )AP. The next step is to investi- 
gate the phase boundary of the TTT sample with respect 
to modulations associated with r 6 . 

Figure 5 (a) shows J c (I?)-curves of the Tuebingen Tri- 



angle Tiling sample at five different temperatures. There 
are two strong sub-peaks for flux densities B « 0.391? i 
and B & 0.65Bi. Similarly to the triangular lattice this 
substructure for (probably irrational) fractionals of B\ 
repeats itself with a somewhat smaller amplitude be- 
tween B — B\ and 2B\. Also in our previous numeri- 
cal simulations we have found this characteristic double- 
peak structure at approximately the same values as in 
the present experimental study [28]. However, we have 
not found a simple connection between these peaks and 
intuitively understandable vortex configurations. 

Regarding the first matching field, there seems to be 
just the standard single maximum for B = B\ at first 
sight. But a zoom-in as in Fig. 5 (b) indicates a double 
peak structure, where the two peaks are merged into a 
broad one with two kinks at B = B\ and at a slightly 
higher value (marked by arrow) . The vertical dashed line 
in Fig. 5 (b) indicates the position of B — (1 + l/r 6 )Bi 
and matches very well the flux density value of the sec- 
ond kink. The second peak hence can be attributed to the 
flux value, where we have one vortex per antidot and one 
additional vortex per 10/10-circle in the sample. One 
possibility is that the additional vortex is at an inter- 
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FIG. 5: (Color online) (a) Critical current I c vs normalized magnetic flux density B / B\ of a Nb structure with a Tuebingen 
Triangle Tiling lattice of antidots at five different temperatures T; (b) and (c) show I C (B) in more detail at B « B\ and 
2Bi, respectively. Arrows point to the excess flux commensurabilities; dashed and dotted vertical lines mark the flux densities 
B/B 1 = 1 + 1/r 6 ((b), dashed), B/B x = 2 + 1/r 6 ((c), dashed) and B/B x = 2 + 2/r 6 ((c), dotted). 



stitial position somewhere inside the circle or that the 
additional vortex pushes the vortex out of the center an- 
tidot, so that the two of them symmetrically occupy the 
circle interior as interstitials. However, we favour a third 
possibility. 

The idea is connected to a recent imaging experiment 
of vortices in a quasiperiodic Penrose pinning lattice [20] , 
which revealed a symmetry-induced formation of giant 
vortices. In the Tuebingen Triangle Tiling situation we 
think that the antidot in the center of the 10/10-circles 
is occupied by a double- vortex at B = (1 + l/r 6 )Bi, 
while all other pinning sites are occupied by one vortex. 
This double-vortex is supposed to compensate the mag- 
netic pressure from the surrounding vortices on the ten 
vortices pinned at the circle line. It stabilizes the whole 
configuration this way. Exactly this vortex configuration 
with ten vortices forming a circle and one double- vortex 



in its center was found in Ref. [20] for pinning sites ar- 
ranged in the Penrose tiling. It is remarkable that in 
the Penrose pattern the pinning site circle is not even 
complete, but has three pinning site free positions. The 
fact that the giant vortex configuration forms anyway, 
strongly indicates its energetic favorability. This should 
be even more favorable in the Tuebingen Triangle Tiling. 
Figure 6 (a) depicts the Tuebingen Triangle Tiling and 
possible vortex configurations corresponding to B — B\ 
(Fig. 6 (b)) and B = (1 + l/r 6 )B 1 (Fig. 6 (c)). 

Similarly to the phase boundary of the Shield Tiling, 
we also find interesting features in I C (B) of the Tuebingen 
Triangle Tiling sample around B = 2£?i, cf. Fig. 5 (c). 
We note the surprising absence of a maximum at B = 
2Bi. Nevertheless we find a clear double peak structure 
similar to the one around the first matching field. The 
first of the two peaks is approximately at a value B = 
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FIG. 6: (Color online) (a) Sketch of a part of the Tuebingen 
Triangle Tiling with colored 10/10-circles and 9/10-circles. 
(b)-(e) Possible vortex configurations related to the critical 
current maxima at (b) B — B±, (c) B — (1 + r 6 )Bi, (d) 



B = (2 + r b )Bi and 
order giant vortices. 



[e) B > (2 + r b )B! including higher 



(2 + l/r 6 )Bi (dashed vertical line) and the second at 
a somewhat higher value, but not exactly at B = (2 + 
2/r 6 )Bi (dotted line). The interpretation of the first 
peak in terms of vortex configurations is obvious. Here 
we have the situation that all antidots are occupied by 
two vortices, whereas those antidots in the center of a 
10/10-circle are occupied by three vortices, cf. Fig. 6 (d). 
The absence of the 2Bi peak indicates that a third vortex 
in the center antidot might be necessary to stabilize the 
configuration. In contrast to the situation around B = 
B\ , it is likely that some vortices from the circle lines sit 
in the interior as interstitials rather than being pinned 
at an antidot. We think this situation is similar to that 
in the Shield tiling, where the caging peak gets stronger 
than the 2B\ peak with decreasing temperature. 

Analogous to the peak at B = (1 + l/r 6 )Bi we should 



find a maximum for the situation when the antidots in 
the 10/10-circle center are occupied by twice the num- 
ber of vortices as the surrounding pinning sites, i.e. by a 
quadruple-vortex. But this seems not to be the case, as 
the second kink of the double peak is found at a higher 
value than B = (2 + 2/r 6 )B 1 (dotted vertical line). We 
believe, that in this second peak another structure in 
the Tuebingen Triangle Tiling comes into play, which we 
call the 9/10-circle. This structure is almost identical 
to the 10/10-circle but with one of the ten circle anti- 
dots sitting at a closer position to the one in the center. 
In Fig. 6 (a) three 9/10-circles are visible. They are al- 
ways directly neighboring a 10/10-circle forming an 8- like 
structure. One possibility to get a stable vortex configu- 
ration is sketched in Fig. 6 (e), where the center antidots 
of the 10/10-circles are occupied by quadruple- vortices, 
the center antidots of the 9/10-circles are occupied by 
triple-vortices and all other by a double-vortex. Unfor- 
tunately we cannot give an analytical expression for the 
corresponding flux density value, as we were not able to 
determine the ratio of the number of 9/10-circles over the 
total number of antidots. 

Probably one has to consider more complicated con- 
figurations including interstitial vortices and also present 
8/10-circles to find the true vortex configurations. Fur- 
thermore, there seem to be some more small shoulders 
and peaks above 2Bi, which we cannot unambiguously 
attribute to simple vortex configurations. 



D. The Transition Temperatures 

We finally consider selected horizontal phase boundary 
T C (B) slices of the three samples. We define T c by a fixed 
ratio R c /R n , where R n is the normal state resistance 
slightly above the superconducting transition at B = 
and R c is a resistance value within the transition regime. 

The transition temperatures of the three samples are 
comparable, cf. Fig. 7 (a). However, the absolute values 
as well as their ratios depend on the resistance criterion. 
For instance, the sample with the triangular antidot lat- 
tice shows a clear two-step (or even three-step) transition 
R(T) for an unknown reason and depending on the crite- 
rion has a smaller or larger critical temperature than the 
ST sample. Also, the general shape of the transition (e.g. 
the curvature of the curve at R c /R n = 0.5) slightly varies 
from sample to sample and all transitions show smaller 
or larger steps. Fingerprints of these stepwise transitions 
can also be found in the T C (B) curves, which are shown 
in Fig. 7 (b)-(d). 

Concerning our discussion of excess flux stabilities, the 
most important observation in the T C (B) data is that 
these commensurabilities for the Shield and the Tuebin- 
gen Triangle Tiling, which we have already found in the 
I C (B) patterns, are also clearly visible mT c (B). All com- 
mensurability features are rather weak for the highest 
resistance criterion R c /R n = 0.7, but get stronger and 
more pronounced with decreasing R c . One exception 
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FIG. 7: (Color online) (a) Resistance R vs temperature T 
close to the transition temperature T c of the three inves- 
tigated samples; (b)-(d) Transition temperature T c vs nor- 
malized magnetic flux density B / B\ for four different resis- 
tance criteria (from top to bottom in each graph R c /R n — 
0.7,0.5,0.3,0.1) of the three investigated samples; R n was 
determined at T — 8.6 K. 



with yet unclear origin has been observed for the Shield 
Tiling, where the first matching peak seems to get weaker 
again for the lowest criterion R c /R n = 0.1. When R c 
is decreased further (and with it the temperature), this 
peak gets stronger again, cf. also I C (B) in Fig. 4, where 
the temperatures are lower than in Fig. 7 (c). This be- 



haviour of the first ST matching peak strength might 
be connected to the transition from a superconducting 
wire network to a thin film with holes, but the observed 
variations and steps in the transition as well as the na- 
ture of the problem confront us with several complica- 
tions for a detailed interpretation and quantitative anal- 
ysis. Due to a not well-defined transition temperature, 
we can not reliably calculate the values for the charac- 
ter ictic length scales such as penetration depth A(T) and 
coherence length £(T), which diverge at T c and hence are 
strongly dependent on the measurement temperature as 
well as on the transition temperature. It is thus not pos- 
sible to reliably calculate, at which temperature the co- 
herence length £(T) is comparable to the width w of the 
superconducting material in between the holes, i.e. at 
which temperature the transition from a superconduct- 
ing wire network (£ > w) to a film with holes (£ < w) 
takes place. For the case of a quasiperiodic arrangement 
of antidots, this transition is not even well-defined, as w 
is spatially varying. 

Nevertheless, for having an idea of the relevant num- 
bers, we can calculate the approximate temperature be- 
low T c , at which the distance between the holes is compa- 
rable to the coherence length. Assuming T c = 8.5 K and 
w ~ 500 nm we find with the expressions from Sec. II 
that f > 500 nm if T > 8.494 K. In other words, as long 
as we are farther away from T c than 6mK the samples 
can be viewed as superconducting films with holes rather 
than as wire networks. 

Although the numbers here are only rough estimates, 
all I C (B) curves above can be considered to have been 
taken at temperatures well below T c and so the observa- 
tions correspond to Abrikosov vortex physics. However, 
due to the differences in the shape of the R(T) transition 
and the absolute values of T c we cannot compare the per- 
formance of different pinning lattices concerning critical 
currents, neither at absolute values of T nor at reduced 
values T/T c . 



V. CONCLUSIONS 

In conclusion, we have experimentally investigated the 
I C (B) — T C (B) phase boundary of niobium thin films 
with quasiperiodic microtopology. We have done this by 
means of transport characterization measurements and 
discussed the results with the focus on observed com- 
mensurability effects for flux densities larger than the 
matching flux. Our experiments confirm some theoreti- 
cally predicted and yet experimentally unobserved excess 
flux matching effects for antidots arranged correspond- 
ing to the Shield Tiling and to the Tuebingen Triangle 
Tiling. For the Shield Tiling, we attribute this excess flux 
matching to an effective caging of Abrikosov vortices in 
quasiperiodically distributed antidot-free hexagons. We 
find similar excess flux commensurability features for an- 
tidots arranged in the Tuebingen Triangle Tiling. How- 
ever, apart from the possibility of caging, the results 
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can be understood as a consequence of recently reported 
symmetry-induced formation of higher order giant vor- 
tices. In particular, we find in both cases that the ad- 
ditional vortices, which create the excess flux commen- 
surabilities, are needed to maximize the local stability 
of the vortex lattice. To finally confirm our interpreta- 
tions of the matching effects and to decide which vortex 
configurations are correct, one needs to carry out imag- 
ing experiments such as magneto-optical imaging, Bitter 
decoration or Hall probe microscopy. 
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